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Abstract 


One-point and two-point exponential functions have been developed and proved to be very effective 
approximations of structural response. The exponential has been compared to the linear, reciprocal and 
quadratic fit methods. Four test problems in structural analysis have been selected. The use of such 
approximations is attractive in structural opumization to reduce the numbers of exact analyses which 
involve computationally expensive finite element analysis. 


1. INTRODUCTION 

The use of detailed, computationally expensive, finite element models has motivated researchers to 
develop approximations of structural response. These approximations are useful for re-design particularly 
with use of optimization techniques, where the number of finite element analyses can be significantly 
reduced. The problem considered here is to construct local approximations using function values and 
derivatives of the structural response at one or two design points. The term local approximation used here 
means that the approximation is valid only in the vicinity of the current design point and is different from 
global approximation methods based on simplified design models or reduced basis techniques which seek to 
approximate the response in the entire design space. 

Specifically, let x be the current design point, where x = (xi, X 2 , ... , x n ) T is a design variable 
vector. Let g (x) be a structural response such as element stress or fundamental frequency, which enters as 
a constraint function in an optimal design formulation. The problem is to construct a local approximation, 
ga(x), based on the function value and derivatives evaluated at x° and possibly another design point. Then, 
subsequent evaluations of the structural response in the neighborhood of x® can be estimated using g a rather 
than the exact response g which will involve finite element computations. A variation of this problem is as 

follows: Let p be a direction vector in design space which has been determined to be desirable in terms of 
reducing the cost function subject to constraints. Usually, p is determined by solving a linear program or 
quadratic program in optimization algorithms. Now, let x^ be a second design point along p such that 

||x — represents a move limit along p. The problem is now to develop a (local) line approximation g a (x) 
such that g a (x) * g(x) for points x along the line joining x° and x 1 , given by 

x = (1-0 x° + C x 1 , 0 < £ < 1 (1) 

Here, the approximation g a is to be constructed using structural response information at x° and possibly x 1 . 

A comparison of various approximation methods has been carried out by Haftka, et. al. (ref. 1) and 
Haftka (ref. 2). The methods include linear and quadratic Taylor series expansions involving first order and 
second order sensitivity analysis (refs. 3-5), approximations based on use of reciprocal design variables 
(refs. 6, 7) and convex approximations (ref. 8). Recently, force approximations have been used by 
Vanderplaats (ref. 9). Use of rational polynomials may be found in Ref. 10. In this paper, exponential 
approximations of the form 


g a = cn x a ‘ 

i=1 ‘ (2) 

are considered and compared with linear, reciprocal and quadratic polynomial methods. It is noted the 
exponential approximation discussed in Ref. (1) is of a different form than that in (2). The motivation for 
choosing exponential approximations of the form in (2) is discussed below. 
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2. BASIS FOR EXPONENTIAL APPROXIMATIONS 


The motivation for approximating structural response using the exponential form in (2) is discussed 
in this section, as also the basis for use of reciprocal variables and force approximations. The basis for 
most approximations comes from the equation 


<*( a ) = t 

A 


(3) 


which states that stress = element force/area. Area A is the design variable here. 

Reciprocal Variables 

The choice of reciprocal design variables is natural, since choosing x = 1/A as a variable results in 
a = c(1 / x) being linear in x: 


g(x) = P x 


(4) 


In the x-space, larger more limits can be imposed on changes in design, leading to faster convergence. 
Now, in a statically indeterminate truss, the stress function is of the form 


o(A) = 


P(A) 

A 


(5) 


The force P is no longer a constant, but dependent on design. The choice of x=l/A is still beneficial as it 
tends to linearize the stress function. In general, a first order Taylor service expansion of g(x) in the 
reciprocals of the variables yi, = 1/xj, i=l, ... , n, written in terms of the original variables, xj, is given by 


( V II 

g a (x) = g(x°) + I (xj -x io ) 
i=1 


MO 


V x i J 


dg/dxj 


(6) 


Force Approximations 

The idea here is to approximate P(A) in (5) by Taylor series as opposed to g(A), and obtain 


ct a (A) = 


P(A 0 ) + -^(A 0 ) (A — A 0 ) 
cLA 


(7) 


In the case when P is a constant, the approximation yields G a = P/A which is exact. Otherwise, curvature 
information is retained in (7) and yields a superior approximation to the conventional tangent approximation 
G a = g(A q ) + da/d A • (A-A 0 ). 

Exponential Approximations 

The approximation introduced in this paper is now discussed. Equation (3) may be re-written as 

G (A) = P A" 1 (8) 
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Thus, the stress is seen to be exponentially related to design variable A. This is the basis for approximating 
structural response in n-dimensional space by 


g a (x) = cflxf i 

i=1 

where xi, X 2 , .... x n are non-negative design variables. Choice of constants C and aj are discussed in the 
next section. 


A second and more general basis for exponential approximations lies in the concept of 'elasticity', a 
quantity used by economists and also relevant in nonlinear stress-strain constitutive laws. Consider a 
function g = g(x), where x > 0 is a scalar variable. The elasticity of the function is defined as 


_ d (In g) 
g d (In x) 


or. 


- _dg/g 
g_ dx/x 


(9a) 


(9b) 


Physically, elasticity may be considered to be in the limit, the percentage change in the function due to a 
percentage change in the variable. For instance, g=x 3 has a value e g =3, and g=px' ] has e g = -1. The 
exponents aj in (2) may be considered to be estimates of the elasticity at the current design point. 

In this section, reciprocal and force approximation methods have been introduced using the 
fundamental equation a = P/A as a basis. Work is being done to generalize these methods to be applicable 
to frames and certain elasticity problems as well. The exponential method of approximation has both o=P/A 
as a basis as well as the concept of elasticity of a function. One advantage of exponential approximations of 
the form in (2) is that, for C > 0, the function ga is a monomial, which opens up the possibility of geometric 
programming (ref. 1 1). 


3. CONSTRUCTION OF THE EXPONENTIAL APPROXIMATION 

The problem is to find the coefficient C and exponents aj, i=l, ..., k, such that the approximate 

n 

function g a (x) = C n*i 1 closely matches the exact function gQO in a neighborhood of the current point x°. 
i=1 

One-point and two-point approximations will now be given. 

LaPoint Approximation 

Here, constants C and {a;} are determined using information only at one point x®. The technique is 
based on matching the function value and shapes of g a and g. This technique has been used in the context 
of unconstrained geometric programming where general functions are reduced to posynomial form. Morris 
(ref. 1 1) discusses an application of this concept to structural design problems. We have, upon takins 
logarithms. 
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( 10 ) 


n 

In g a = In C + £ aj In Xj 
i=1 


Differentiating with respect to xj yields 


3g a /3xj = g a 



( 11 ) 


Note that ga(x^) = g(x®). Equating 3ga/3xj in (1 1) to the exact slope 3g/3xj at x® yields the exponents 



dg / 8xj 
g 


( 12 ) 


The coefficient C is then obtained from ga(x®) = g(x®) as 

c=g/nx?- |x G 

i 


(13) 


2-Point Approximation 

Information at two points are used to construct the exponential approximation. Let x° be the current 
design point and x 1 be a second point, which usually is a point along a desired search direction in design 
space. The quantities g(x°), Vg (x°) and g (x 1 ) are now used to determine C and {aj}. A least squares 
formulation is adopted herein. The variable C and {aj} are obtained from the minimization problem 


E 


( 


n 


A 


g 0 -c nx i( j 

i=1 J 


gi-cnxJi 

, i=1 


+ z 

j=1 


^gp 

a Xj 


\2l 


— C n x jo • In x 


i=1 




(14) 


The minimization of the least squares objective function E is carried out using a modified Newton algorithm, 
with a Levenberg-Marquardt correction to the Hessian when descent is not obtained (ref. 12). The 
algorithm requires the gradient vector 


VE = (8E / 3C, 8E / 8a-|, ••• , 8E /8a n ) 


and the Hessian 


a 2 E/ac 2 


h e = 


symmetric 


8 2 e / acaa-j 

a 2 E / 3a 2 ••• 


a 2 E/acaa n 

3 2 E / 3ai 3a n 
3 2 E / a 2 


(15) 


(16) 
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These derivatives are computed from analytically derived expressions. The least squares algorithm is given 
below. 


Algorithm 2-Point Exponential 

Step 1. Choose the initial estimates of C and {aj} from (12), (13), and e Q = 0.001 
Step 2. Solve 

(H E + 8l)5 = -VE T 


(17) 


and update Cnew ~ C + 8i, (aj)new — aj + 5j+i, i — 1, ... , n. 


Step 3. Evaluate E nC w If E n ew < E the set C = C n ew» {aj} = {aj) new , reduce e, say, e = e/ 10 (if e < 
E 0 , set E = £ 0 ) and go to step 2. If E new > E, then increase £ = 10.E and go to step 2. 

The procedure above is terminated when relative and absolute reductions in E for three consecutive iterations 
are less than a specified tolerance. 


4. TEST PROBLEMS AND RESULTS 

Four test problems relating to structural design have been considered. The 1 -point and 2-point 
exponential approximations developed in Section 3 are examined. Comparison of the approximation to the 
original function is done along a line joining two design points x°, x 1 , or at points x where 

2Si = 0 — C) 2^ + C 2^ 

where C, is scalar variable, 0 < £ < 1 . For comparison, the linear (tangent) approximation based on 

g a (x) = g(x°) + Vg(x°)*(x-x°) ( 15 

the reciprocal-linear approximation given in (6), and the quadratic polynomial along the line given by 


ga(0 = a + bC + cC 2 (2C 

where coefficients a, b, c are obtained from g(x°), g(x ] ) and dg/d£ (at £=0) = Vg(x°) • (x L x°). Thus, the 
1-point exponential, linear, and reciprocal require only g(x°), Vg(x°), while the 2-point exponential and 
quadratic polynomial require, in addition, g(x J ). The error between g and ga along the line is shown both 
graphically as well as quantitatively through a relative error criterion 


RELEER = 



[(gi ~ gal) / giT 


and a maximum error criterion 


( 21 ) 
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MAXERR = max |g, - g^l 
1 < i < N 


(22) 


Above, gi = g(x (£j)) is the exact function evaluated at the i* discretization point along the line in (18), g a i is 
the approximate function evaluated at £i, and N, the number of discretization points, is chosen equal to 20. 

Cantilever Beam 

The axial stress function in a cantilever beam of rectangular cross section, subjected to axial and 
transverse loads, is given as 


, , 1000 6000 
o(x) = + ■ 


xix 2 


V Y 2 
X 1 X 2 


(23) 


where xi, X 2 are the width and depth of the cross section, respectively. The choice of design points is x° 
= (1,2 ) t in., x 1 = (5,8) T in. 

Referring to Fig. 1, the 2-point exponential is in excellent agreement with the original function. The 
1-point exponential behaves just as well as the 2-point exponential and is not shown in the figure. The 
exponential approximation to o(x) in (23) is of the form 


o a (x) = 6727.2 xf°' 891 xjj 1 ' 750 


(24) 


The quadratic polynomial (Fig. 1), as well as the tangent and reciprocal approximations behave very poorly. 
The values of RELEER and MAXERR in (21), (22) for this problem are given in Table 1 . It is noted that 

various choices of x° and x 1 have shown the same trend. 


Table 1. Cantilever Beam 


Approximation Method 

Relative Error 

Maximum Error 

1-ooint exnonential 

0.493 

13.1 

2-point exponential 

0.492 

13.0 

Linear (taneentl 

700.0 

0.165 E5 

Reciprocal 

105.4 

0.022 E5 

Quadratic Polvnomial 

96.9 

0.033 E5 
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Tension-Compressio n Spring 

The shear stress function in a spring design problem, with X! = coil diameter and x 2 = wire diameter 
is given by 


T(x) = 


8000 x-| 
7C X 2 


4xi - X 2 

v 4 Xi -4x 2 ; 


+ 0.615 — 
X 1 


(25) 


For this problem, x° - (1.0, 0.3) T and x 1 = (0.3, 0.05) T in. As with the beam, the 2-point and 1 -point 
exponentials are also in close agreement with the original function. Table 2 provides an error summary for 
all the methods. The linear, reciprocal and quadratic polynomial are poor by comparison (Fig. 2). Other 
choices of x°, x 1 show the same trend for this problem. 


Table 2. Tension -Compression Spring 


Approximation Method 

Relative Error 

Maximum Error (x \(P) 

1 -point exponential 

0.091 

0.496 

2-point exponential 

0.088 

0.067 

Linear (tangent) 

2.464 

7.265 

Reciprocal 

1.640 

5.876 

Quadratic Polvnomial 

11.060 

3.539 


Three Bar Symmetri cal Truss 

The natural frequency of a three bar truss (ref. 3) with x j, x 2 = cross sectional areas, is described 
by the function 


(o(x) = — 

2 V 2 Xi + x 2 

Two sets of design points, leading to different performances, are chosen. These sets are 
I- x° = (3,4)Ti n 2 > x 1 =O0,5)Ti n 2 

II. x0 = (5,5)Ti n 2 , x 1 = (1,10) T in 2 


(26) 


(27) 


SeLl: Referring to Fig. 3 and Table 3, the 2-point exponential based on the best fit formulation yields the 
best approximation, with 


0.261 „ -0.342 


co(x) = 0.2877 x p D1 x 2 


(28) 
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The 1 -point exponential is poorer, with 


<o(x) = 0.2635 x?- 32 x 2 °- 32 ( 29) 

Set II : Along the search direction, the original function is quite flat. In fact, the 1-point exponential 
provides a relatively poor approximation because of the flat nature of the function. The quadratic polynomial 
is best here. Even though the 2-point exponential is second-best, (Fig. 4), the best- fit nature of the 
approximation, while averaging the error, does not provide an interval within the line where the error is 
small. This may cause difficulty for designs near the optimum. Finally, the use of reciprocal variables does 
not show any advantage over the direct variables for this case. 


Table 3. Three Bar Symmetrical Truss 


Approximation Method 

(Set 

Relative Error 

I) 

Maximum Error 

(Se 

Relative Error 

til) 

Maximum Error 

1 -point exponential 

0.247 

0.0286 

1.386 

0.065 

2-point exponential 

0.059 

0.0053 

0.605 

0.028 

Linear (tangent) 

0.810 

0.100 

1.210 

0.060 

Reciprocal 

0.199 

0.022 

1.902 

0.060 

Quadratic Polynomial 

0.128 

0.012 

0.074 

0.004 


Ten Bar Truss 

The ten cross sectional areas of the truss shown in Fig. 5 are the design variables. Points x° and x 1 
are chosen as the initial and optimum design obtained in Ref. (1), as 

x° = (5., 5., 5., 5., 5., 5., 5., 5., 5., 5.) T in 2 

and (30) 

x 1 = (7.94, 0.1, 8.06, 3.94, 0.1, 0.1, 5.74, 5.57, 5.57, 0.1) T in 2 

Again, both 1 -point and 2-point exponentials provide excellent approximations. The number of 
design variables do not seem to affect their quality. Interestingly, the reciprocal approximation provides 
equally good results, but to within a certain distance from x®. Near x 1 , the reciprocal abruptly diverges 

(Fig. 6). With smaller move limits, of course, the reciprocal will be excellent for this problem. The 

quadratic polynomial provides a good approximation for this problem. Error estimates are given in Table 4. 
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Table 4. Ten Bar Truss 


Approximation Method 

Relative Error 

Maximum Error 

1 -point exponential 

0.088 

1.528 

2-point exponential 

0.021 

0.204 

Linear ftanaentl 

0.689 

7.907 

Reciprocal 

0.599 

13.980 

Quadratic Polynomial 

0.081 

0.803 


5. CONCLUSIONS 


Exponential functions of the form CflXj' have been used to approximate structural response. Both 

1-point and 2-point approximations have been used to determine C and { ai }. The 1 -point involves matching 
function and derivative values at the current design. The 2-point method is based on minimizing a least 
squares function by modified Newton's method. The basis for exponential approximations is from two 
sources: one is from structural theory, where o = P/A can be written as a = PA -1 , while the other is from 
economics, where a function g = cx a has an elasticity equal to the exponent a. The restriction of exponential 
approximations is xj > 0. An advantage is that the approximating function is valid for any type of structure 
or type of structural response. Further, the exponential approximations when applied to the cost and 
constraints of an optimal design problem have the potential for being used in conjunction with geometric 
programming which can effectively solve the subproblem. 


Results on three out of the four structural problems considered have shown that the exponential 
functions have provided excellent approximations, with essentially no error, even for large distances in the 
design space. The linear, linear-reciprocal and quadratic polynomial are much inferior. 

On one of the problems involving natural frequency of a 3-bar truss, it is observed that the function 
is essentially flat or linear. In this case, the 2-point exponential (based on a best-fit) does better than the 1- 
point, but the quadratic polynomial is superior, Thus, for linear or nearly linear functions, the linear 
approximation is to be preferred. For other cases, the exponential has shown to be a powerful method of 
approximation. 
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Figure 5. Ten Bar Truss Problem. 
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